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Abstract 

In grand unified theories with gauge groups larger than SU(5), the 
multiplets that contain the known quarks and leptons also contain 
fermions that are singlets under the Standard Model gauge group. 
Some of these could be the dark matter of the universe. Grand unified 
theories can also have accidental U(l) global symmetries (analogous 
to B — L in minimal SU(5)) that can stabilize dark matter. These 
ideas are illustrated in an SU (6) model. 

1 Introduction 

It seems a strange coincidence that the cosmic densities of dark matter and 
ordinary baryonic matter are of the same order of magnitude, given that in 
most theoretical scenarios they are generated by unrelated mechanisms in- 
volving different particles, forces, and parameters. This coincidence suggests 
that the dark matter and baryonic matter may have been "cogenerated" in 
the early universe, i.e. that the dark matter is a product of the same pro- 
cesses that created the cosmic baryon asymmetry. There is a rapidly growing 
literature studying various ways that this might have happened [1, 2, 3, 4, 5]. 
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The first papers to propose this possibility [1] were based on the idea that 
primordial asymmetries in baryon and lepton number (B, L) were partially 
converted into an asymmetry in some other global quantum number (call 
it X) by sphaleron processes [6] when the temperature of the universe was 
above the weak interaction scale Mw Assuming X to be conserved (or nearly 
so) at low temperatures, the lightest particles carrying this quantum number 
would be stable and could play the role of dark matter. What would result 
from such a scenario is "asymmetric dark matter" [7]. Many other scenarios 
for generating asymmetric dark matter have been proposed [2, 3, 4, 5]. In 
some of these scenarios ordinary matter and dark matter are converted into 
each other by perturbative processes involving higher-dimension operators 
[2]; and in others by sphalerons (or by both sphalerons and higher-dimension 
operators) [3]. In some scenarios, the dark matter carries baryon number 
which compensates for the non-zero baryon asymmetry of ordinary matter 
[4]. And many papers propose still other mechanisms [5]. 

What we suggest here is the possibility that not only are dark matter and 
ordinary matter "cogenerated" but that they are "unified" in the sense of 
grand unification. The point is that grand unification naturally supplies sev- 
eral of the ingredients needed for the generation of asymmetric dark matter. 
First, grand unification based on groups larger than SU(5) involves fermion 
multiplets that contain non-Standard-Model fermions that could play the role 
of dark matter. In particular, in SU(N) with N > 5, the quark-lepton multi- 
plets contain several fields that are singlets under the Standard Model group 
G SM = SU(2) L x SU(3) C x U(1) Y . In SU(6) or SU(7) models with three 
families of quarks and leptons, for example, anomaly-free sets of fermion mul- 
tiplets must contain at least six Standard-Model-singlet fermions. For larger 
groups the number grows rapidly. 

Second, it is not uncommon in simple unified models for there to be global 
quantum numbers that are "accidentally" (or "automatically") conserved, 
just as B — L is accidentally conserved in the simplest SU(5) model. Such 
quantum numbers could play the role of X that stabilizes the dark matter 
particles. We shall illustrate these ideas in a simple SU (6) x Zn grand unified 
theory (GUT). 

Another interesting feature of large unified groups is that they can contain 
additional non-abelian factors at low energies (besides those of Gsm), whose 
sphalerons could convert baryons and leptons into dark matter particles; but 
we shall not explore that possibility in this paper. The illustrative model 
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that we shall describe uses perturbative processes to convert matter and 
dark matter into each other [2]. 

2 An SU(6) Model 

We shall now present the details of the model. Its symmetry group is SU(6) x 
Zn, where N may be any integer greater than 4. The fermions of each quark- 
lepton family consist of the anomaly-free set of SU(6) multiplets 15 + 6 + 6 
plus three SU (6) singlets. These are shown in the left columns of Table I. 
The fundamental indices of SU(6) are denoted by the capital latin letters, 
A, B, C, etc, which run from 1, 6. Here and throughout the paper, we 
suppress family indices. For those fermion multiplets that contain the known 
quarks and leptons we use capital The fermion multiplets denoted by the 
letters ip and ( and rj contain only new fields. 



Table I: The fermions of a family, and the Higgs fields. 
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The right columns of Table I show the Higgs fields and their SU (6) x Z N 
transformation properties. The known fields all transform trivially under 
Zn- The gauge symmetry breaking occurs in three stages: 

(1) At the unification scale M GUT , SU{6) is broken to [SU(3) C x SU(2) L x 
U(1) Y ] x U(l) 6 (which is contained in the SU(5) x Z7(l) 6 subgroup of SU(6)). 
This breaking is done by an adjoint Higgs field fig, whose vacuum expecta- 
tion value (VEV) points in a direction that is a linear combination of the weak 
hypercharge generator Y/2 = diag(^, |, — |, — |, — |, 0) and the U(l)e gener- 
ator T 6 = diag{— |, — -|, — |, — |, — |, 1). We shall always denote an SU(5) 
index (which takes values 1,2,3,4,5) by a, (3, etc.; an SU(2)l index (which 
takes values 1,2) by i, j, etc.; and an SU(3) C index (which takes values 3,4,5) 
by a, b, etc. Thus, for example, ip A = (^ a ,^ 6 ) = i>a, V'e)- 
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(2) The U(1)q is broken at a scale M', which is somewhat larger than a 
TeV, by the fundamental Higgs multiplet Ha, whose VEV points in the 6 
direction, i.e. (if 6 ) ~ M'. We assume that the mass of H 6 is of order M', 
but that its other components all have mass of order Mqut- (This is the 
usual kind of split-multiplet problem of GUTs, analogous to the well-known 
"doublet-triplet splitting problem".) Below M', the gauge group is just the 
Standard Model group Gsm = SU (3) c x SU (2)lxU (l)y . There is, of course, 
a Z 1 gauge boson whose mass is of order M' that couples to T 6 . 

(3) The electroweak breaking is done by the two Higgs multiplets denoted 
by 0, which contain the SU(2)l doublets <pi and 0* 6 . We assume that one 
linear combination of these two doublets is tuned to be light (i.e. of order the 
weak scale in mass) and is the Standard Model Higgs doublet $j that obtains 
a VEV, while all the other components of 4>a and $ AB } have superheavy mass. 
(In particular, the colored components, which can mediate proton decay, have 
mass of order M GUT .) 

Under the subgroup £77(5) C £77(6), the non-singlet fermions of Table I 
decompose as follows 
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The Standard Model quarks and leptons are the \l>[ a/3 l = 10 and \I/ a = 5. 
The extra 5 + 5 pair of £77(5) will "mate" to get O(M') masses. Specifically, 
the 5 = ^ [a6] will obtain mass with 5 = ip' a through a term (^ a ^ip' a ){H 6 ), 
as will be seen. There are also two £C/(5)-singlet (and thus Gsm- singlet) 
fermions in the multiplets shown in Eq. (1). In order for these to get mass, 
we introduce gauge-singlets fermions denoted r\ and ( that will mate with 
them to get Dirac masses. 

The most general renormalizable Yukawa couplings allowed by SU (6) x 
Z N are of the following forms (we suppress family indices): 
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C'Yukawa — &SM ~\~ £ pp ~\~ £ singlet 

C SM = Y u (^ AB ^f CD ) ( f> EF ^+Y d (¥ AB ^ A ) ( j )B + Y^a^^ + Mr^) 
C fT = h{¥ AB H' A )H B 

Csinglet = / 2 (^M)#* A + hW A QH* A + WaVW* 

(2) 

The first three terms in C SM have the effect of coupling the Standard Model 
Higgs doublet (which is a mixture of <pi and 0* 6 ) to the known quarks and 
leptons (which, as noted above, are contained in the multiplets that are 
denoted by capital \l/'s). The term J~>p~p = fi(^ AB ty' A )HB gives O(M') mass 
to the "extra" 5 + 5 pair of fermions, as already mentioned. The terms in 
£ S ingiet couple the Standard- Model-singlet fermions \l/ 6 and tp^ (see Eq. (1)) 
to the gauge-singlet fermions denoted rj and ( so that all the singlet fermions 
can get mass. All these Yukawa terms and the masses that come from them 
will be examined in more detail shortly. 

The most general Yukawa terms allowed by S77(6) x (shown in Eq. 
(2)) and the most general Higgs potential allowed by 577(6) x Z N (which, 
incidentally, includes terms such as 4> AB (f) A Hc^lB) happen "accidentally" to 
be invariant under a global U(l) symmetry, whose generator we will call T. 
The T charges of the various multiplets in the model follows (given 

in parentheses): ^(1), ^ A {-\), tf(0), ^(-|), ^(3), C(6), <p [AB] {-2), 
4>a{-\), H a {\), and H|(0). This U(1) T symmetry is unbroken by GUT- 
scale VEVs, but is spontaneously broken at the scale M' by (H 6 ), which also 
breaks the gauge group C/(l) 6 , leaving an unbroken global symmetry U(l)x, 
whose generator is given by 

X = \T + |T 6 . (3) 

This generator X will play a crucial role in what follows as the quantum 
number that stabilizes dark matter. It is analogous to the conserved global 
symmetry B — L in minimal SU (5) models, where B — L is a linear combina- 
tion of a global charge that is accidentally conserved by the Yukawa couplings 
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and a gauge generator (specifically, the weak hypercharge Y). Here, both X 
and B — L are conserved by low-energy couplings and VEVs. 

All the Higgs fields in Table I have 0(M GUT ) masses except (a) the Stan- 
dard Model doublet $j, which is a linear combination of the doublets fa and 
0* i6 ], and has mass of order 100 GeV, and (b) the Standard- Model-singlet 
Higgs Hq, which has mass of order M' > 1 TeV. We will call $j and H$ 
"light Higgs fields" to distinguish them from "superheavy Higgs fields". It 
is easy to check explicitly that the light Higgs fields $j and H e are neutral 
under X (as must be the case, of course, if X is left unbroken by their VEVs). 
Therefore, when fermions absorb or emit these light Higgs fields they do not 
change their X charge. The same is true of the emission and absorption of 
light gauge bosons (i.e. those of SU(3) C x SU(2) L x U(1) Y x U(l) 6 ). 

In other words, except through extremely slow processes mediated by 
bosons with superheavy masses, fermions do not change their values of X. 
It is therefore very useful to classify the fermions of this model by their X 
charges. This is done in Table II. 

Note that Table II also lists quantum numbers called B , L , B±, L±, L 2 . 
These are defined as follows: B n = B5\x\ n , L n = L6\x\ n - 111 other words we 
define separate baryon and lepton numbers for each \X\ sector. For example, 
a lepton with X = —2 has L 2 — 1, but Lq = L\ = 0. (It should be 
noted that B = B + B 1 , L = L + L x + L 2 , and X = W X - L 1 - 2L 2 .) The 
reason for defining these baryon and lepton numbers is that the emission and 
absorption of "light Higgs fields" and "light gauge bosons" do not change the 
values of B, L, and X of a fermion and therefore also separately conserve the 
quantum numbers B n and L n . This will be important in our later analysis. 
Later we shall introduce four-fermion operators, generated by the exchange 
of very heavy bosons, that conserve B, L and X but violate B n and L n . Such 
processes will be needed to re-distribute particle asymmetries, i.e. convert a 
primordial asymmetry in one global charge into the other global charges, so 
that matter and dark matter asymmetries will end up being related to each 
other. 
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Table II: The fermions and Higgs fields, classified by their X values (0, 
±1, ±2), where X = |T + ^Tq. The baryon and lepton numbers are defined 
by B n = B5\x\ n , L n = L5\ X \ n - 
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Observe that the fermions with X = are just the known quarks and 
leptons of the Standard model. Since the light Higgs fields have X = 0, their 
Yukawa couplings only couple these Standard Model fermions to each other 
and give them Dirac masses with each other. These couplings come from the 
terms Csm i n Eq. (2). Specifically, the Y u term gives mass to up-type quarks 
via Y u (^ [ab W cl] )((j) [26] ) oc Y u (u c u)v. The Y d term gives mass to down-type 
quarks and charged leptons via Y d (^ a2 ^ a + ^ 2 %i)(0 2 ) oc Y d (dd c + £ + £~)v. 
The Y v term gives the Dirac neutrino masses via Y v (^f-$f) ((f)* 2 ) oc Y v (yN c )v. 
And the Mr term gives the superlarge Majorana masses to the right-handed 
neutrinos: M R (^) = M R (N C N C ). 

As can be seen from Table II, the sector of fermions with X = ±1 con- 
tains (for each family) a 5 and 5's of SU(5), namely and ip' a . These 
"mate" to obtain masses of O(M') via the Yukawa term in C F p, which gives 
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h{^H' a ){H % ). 



The X = ±1 sector also contains (for each family) the singlet fermions 
^ 6 with X = -1, and rj with X = 1. The term f 2 (^ 6 7])H* 6 in £ sin9 / et 
couples them together into massive Dirac particles. (Note that the gauge- 
singlet fermions 77 have been introduced into the model just to give mass to 

which otherwise would remain massless.) The term f^ip^rf)^ in C S i ng iet 
has the effect of mixing these singlet fermions with neutrinos in the 5 + 5, so 
that the neutral fermions in the X = ±1 sector actually have a 2 x 2 mass 
matrix (actually 6 x 6 if one takes into account that there are three families) 
of the following form: 



If all the Yukawa couplings in Eq. (4) were of order 1, then all the masses 
of the fermions in the X = ±1 sector would be of O(M'). There have to 
be, however, some particles to play the role of dark matter. Since in this 
scenario the present number densities of dark matter particles and baryons 
will turn out to be of the same order of magnitude, the mass m DM of the 
dark matter particles should be roughly of order 1 GeV. There are various 
ways this can be the case. One simple way is that the Yukawa couplings that 
we have denoted f 2 in Eqs. (2) and (4) are of order (1 GeV)/M' < 10~ 3 . 
We shall assume this to be true and also assume that the Yukawa couplings 
denoted j\ and f 3 are significantly larger than f 2 . In that case, the lightest 
I/O fermions are the Dirac fermions made up of the gauge singlets ^/q and 
rj. These will be the dark matter particles of the model. We will denote these 
dark matter particles sometimes as (As noted, and as can be seen 

from Eq. (4), these mix with angle 0(v/M') with weak-doublet neutrinos of 

mass O(M'). Thus the dark matter particles have O (jj^r) couplings 

to the Standard Model Z boson.) 

We come, finally, to the X = ±2 sector of fermions. This consists (for 
each family) of a Standard-Model-singlet fermions with X = —2 (namely, 
■0g), and with X = +2 (namely (). The term f^(il)' A C)H* A in C S ingiet couples 
these together to make massive Dirac particles. (Note that the gauge-singlet 
fermions denoted by the letter ( have been introduced into the model just to 
give mass to the X = —2 fermions.) 
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3 The Processes that Redistribute Asymme- 
tries 



As noted before, there must be processes that conserve X and B — L but 
violate B n and L n in order to redistribute asymmetries in quantum numbers 
and thus relate the matter and dark matter asymmetries. Such processes are 
needed for other reasons as well. For example, they are needed to allow the 
colored 1^0 particles (i.e. those with B\ ^ 0) to decay. (Such particles 
if stable and light would have been seen at accelerators, and if stable and 
heavy would contribute too much to the dark matter density of the universe.) 
These ABi ^ decays can be relatively slow, as long as they occur early 
enough not to interfere with primordial nucleosynthesis. 

The processes that we postulate to violate B n and L n are very simple. 
They are given by four-fermion operators of the form (ip' ip' B )(^ ^a), and 
more precisely by 

O x = (M 1 )- 2 (W 6 m^), A(B ,L ,B 1 ,L 1 ,L 2 ) = (0,-1,0,2,-1), 

2 = (M 2 )- 2 (V7 6 ^)(^>e), A(fl , Lo, B u L u L 2 ) = (§, 0, -§, 1, -1), 

3 = (M 3 )- 2 (^V9(^a), A(B ,L ,B 1 ,L 1 ,L 2 ) = (-1,-1,1,1,0). 

(5) 

It is easy to check that these conserve B = B +B 1 , L = L +Li+L 2 and X = 
3Bi — Li — 2L 2 . These operators can arise in a simple way from integrating 
out heavy fields as follows. Suppose there is a boson 0^ and gauge-singlet 
fermions \& and fj. Suppose (f>A, ^, and fj have the same SU (6) x Z^, quantum 
numbers as the particles 4>a, an d f], respectively. We assume that <p a and 
fj have mass of a scale we will call Ma, where Ma ^ M', and that 6 has 
mass of order 1 to 10 GeV, while ^ is massless. (It is also assumed that 4>a 
has vanishing VEV.) With these quantum numbers, these particles have the 
Yukawa couplings 

£a = f(* A *W A + f'^' A vW A - (6) 

Note the similarity to the terms in Eq. (2) with coefficients Y v and f±. We 
can ensure that these are the only Yukawa couplings that <pA, ^, and fj 
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possess, either by assigning them suitable Zn charges, or, even more simply, 

by positing another Zm symmetry under which (pA — > z<Pa, ^ — > z^f, and 

fj — >■ zfj, while all other fields transform trivially. (We can assign the massless 

fermion ^ lepton numbers L — 0, L\ — 1, L 2 — 0, with the light boson c/> 6 

having all baryon and lepton numbers zero.) 

Then the box diagram shown in Fig. 1 gives rise to the operators in Eq. 

(5), with Mi, M 2 , Ms ~ Ma- Note that the operator O3 directly gives the 
2 

decay ip' a — > + ^ + ip' 2 . The initial particle is an X = — 1 antiquark 
with mass of O(M'). The first two final state particles are ordinary Standard 
Model particles, namely an antiquark (d c L ) and a lepton (vl)- The third final 
state particle is i/j' 2 , which mixes with the dark matter particles (^/q,i]) as 
can be seen from Eq. (4). In fact, the operators in Eq. (5) allow all fermions 
of the model with masses of order M' to decay ultimately to ordinary quarks 
and leptons and the dark matter particles (^6,^)- 



^4 * 



Fig. 1 




Figure 1: The diagram that gives the operators 0\ (if A = 6 and B = i), 
2 (if A = 6 and B = a), and £> 3 (if A = a and B = i). 

4 The Cosmological Scenario 

Now let us outline the sequence of events in the early universe that generate 
the current baryon and dark matter abundances in this model. 

Stage 1, which happens when the universe is above a superlarge temper- 
ature Tlg, is the genesis of an asymmetry in lepton number. Often this is 
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assumed to happen through the decays of superheavy right-handed neutri- 
nos (here called \&) [8]. However, as will be seen, this will not lead to the 
generation of any dark matter asymmetry in the present scenario, since the 
right-handed neutrinos have X = 0. We will therefore assume that the pri- 
mordial asymmetry is of fermions that carry both L and X. In particular, we 
shall look at the case in which an asymmetry of the fermions ( is generated. 
How this might happen will be discussed later. 

Stage 2 is the period Ta < T < Ma, during which both SU(2) L sphaleron 
processes and the scattering processes involving the operators Oi, 2 and 
O3 are in equilibrium. The freeze-out temperature Ta of those scattering 

processes is given roughly by Ta ~ Ma ( ^Mp/^ ) ' wnere 9 * s ^ ne effective 
number of massless particle species at Ta- We shall assume that Ma is 
large enough that Ta > M', the scale at which C/(l) 6 breaks and the non- 
Standard-Model quarks and leptons get their mass. During stage 2, the 
initial asymmetry in ( (and thus in L and X) is reshuffled by sphalerons and 
scattering processes among the various particle types, leading to asymmetries 
in all of the quantum numbers Bq, Lq, B\, Li, and L2 that are of the same 
order. These will be computed shortly. 

Stage 3 is the period T dec < T < T A , where T dec is the temperature at 
which the particles with mass of O(M') decay. When the temperature falls 
below Ta, the relative values of B , L , B\, L\, and T 2 freeze. Sphaleron 
processes continue until some temperature T sp h ~ 200 GeV, but do not af- 
fect these ratios, which change again only when the fermions with mass of 

order O(M') decay (out of equilibrium) via the operators 0±, 2 , and C 3 at 

/ \ 1/2 

Tdec ~ M' ( 96 ^2y3 J • K should be noted that when T falls below M' (which 

happens during stage 3, since Td ec < M' < Ta), annihilations (mediated by 
SU(3) C x SU(2) L x U(l)y x £7(1)6 gauge interactions wipe out almost all the 
particles of O(M') except for the asymmetric components. (For example, 
if the asymmetry in B\ is positive, the density of B\ < particles will be 
driven to a value much less than that of B\ > particles.) 

Stage 4 is the period when T < T dec . At this point, the particles with 
mass of O(M') decay out of equilibrium via the operators 1 , 2 , and 3 . 
These decays violate B , L , B 1: T l5 and T 2 , and so again reshuffle the ratios 
of these quantum numbers. In particular, they set Bi and L 2 to zero; but 
they leave a non-zero L\ in the form of the dark matter particles 
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which we will compute below. At this point, the only remaining fermions are 
the known Standard Model quarks and leptons, the dark matter particles 
(\& 6 ,77), and the massless SSf particles (whose contribution to the radiation 
density at the time of nucleosynthesis is equivalent to less than half of a 
neutrino species). 

We assume that Td ec > ™>dm ~ 1 GeV. (This means, for example, that if 
M' ~ 1 TeV, then Ta must be also be of order 1 TeV, while if M' ~ 10 TeV, 
then Ta must be less than about 100 TeV.) When T falls below m DM the dark 
matter particles (^e,r]) start to annihilate with their antiparticles through 
the diagrams shown in Fig. 2. Since the 4>q have been assumed to be light (of 
order 1 to 10 GeV) these annihilations efficiently reduce the number of dark 
matter anti-particles to much below the number of dark matter particles. 
That is, the dark matter is asymmetric. 



^ 6 * 



"06 

^6 ^ 

Fig. 2 

Figure 2: The diagram by which dark matter particles and antiparticles can 
annihilate into massless fermion antifermion pairs: ^6 + ^ — > ^ + 

5 Computing the Relative Asymmetries 

We now calculate the matter and dark matter asymmetries that result in 
this model. During stage 2, all of the fermions of the model except \1/ and 
fj are relativistic, so we can neglect their masses. Call the fermion chemical 
potentials /i qo , fj,£ , /i qi , /i^, and /i^ 2 , for the quarks and leptons having various 
values of X. From the fact that the sphaleron processes and the scattering 
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processes that involve the operators Oi, and O3 are in equilibrium, one has 
the following relations 

sphaleron : = 3fi qo + /^ , 

0\ : = ^-2^+^, (7) 

The equilibrium of the scattering processes involving C 2 does not give an 
independent relation. From Table I, it can be seen that the number of species 
of fermions of each type (i.e. the statistical weight), counting the number of 
families, colors and polarizations, is given by g(qo) = 3 ■ 3 • 4 = 36, g(£o) = 
3-1-3 = 9, g{q{) = 3-3-2 = 18, g(£ 1 ) = 3-1-6 + 1 = 19, and g{£ 2 ) = 3-1-2 = 6. 
(The extra 1 appearing in g{£\) is due to the fermion \&.) In a co moving 
volume of the universe, the asymmetries in the quantum numbers are related 
to the chemical potentials by the relation TVj oc gifiiT 2 R(T) 3 , where R(T) is 
the scale factor of the universe when the temperature is T. Therefore, one 
has 



(8) 



















L 2 





where K depends on the temperature and volume. Therefore, from Eqs. (7) 
and (8), one has that during stage 2 

= \B + |L , 

= + §A) - \B\ - YgLi, (9) 
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We assume that during stage 1 primordial asymmetries were generated in 
the quantum numbers X and B — L. After stage 1, however, these quantum 
numbers are conserved. Thus we have two further relations 

X = 3.Bi - Li - 2L 2 = a, 

(10) 

B — L = Bq + B\ — Lq — L\ — L 2 = 6, 
where a and 6 are constants. Eqs. (9) and (10) can be solved to obtain, 

B » = -9TT9 ( 37a - Qlb ) > 
L = ^(37a-61b), 

B l = 9T19 ( 83fl - 50& ) > (H) 

L 2 = ^Ti9 (-86a + 266) . 

In stage 4, after T has fallen below T deC) the particles with mass of order 
M' decay into the ordinary quarks and leptons of the Standard Model, plus 
the dark matter fields ^q. We will assume for ease of discussion that the 
heaviest O(M') particles are those with Bi ^ followed by those with L 2 ^ 
0. Then the B\ ^ particles will decay via the operator C 3 . From Eq. 
(5) one sees that these decays will change the particle asymmetries in the 
proportions A^o = -ABi, AL = -3A£i, and ALi = 3AB 1 . Therefore, if 
B 1 — > B[ = B 1 + ASi = 0, one has 
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(-86a + 266) . 



The L 2 7^ particles decay via the operator Oi, changing the asymmetries 
in the proportions AL = AL 2 and AL X = —2AL 2 , as can be seen from Eq. 
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B f = 


2 
119 


(a + 86) , 


L f = 




(121a- 1036), 


L lf = 


—a. 





(5). Therefore, if L 2 — > 0, one ends up with the final values of the quantum 
numbers being 



(13) 



It should be noted that the final symmetry in L\ is in the form of the 
massive dark matter particles (^,77), and not in the form of the massless 
fermions (It is easily shown in the following way that there is no asym- 
metry in ^ . One can assign an exactly conserved quantum number Z to \I/ 
and 06; with both these particles having Z — 1 and all other particles having 
Z = 0. By conservation of Z, and the fact that no asymmetry of Z existed 
initially, one has that N$ — = 0. But eventually all the ip G and their 
antiparticles decay by 6 — > ^6 + ^ , which drives N^ 6 and thus Nq, to zero.) 

Let us suppose that the primordial asymmetry generated during stage 1 
is in the number of ( particles. Since these have X = 2 and B — L — 1, it 
follows that a = 26. From this and Eq. (13), one has that 

L lf 119 



(14) 

B of 10 V ; 

This is the present ratio of the number of dark matter particles to the number 
of baryons in the universe. It is rather remarkable feature of models of this 
type [1, 2, 3] that this ratio is predicted, and therefore that the mass of the 
dark matter particle is predicted. Of course, different grand unified models 
would give different predictions. 

Finally, we come to the question of the primordial asymmetry generated 
during stage 1. One could imagine that this asymmetry was in ordinary 
leptons, from the decay of superheavy right-handed neutrinos as has been 
much studied [8]. However, that would give a = and 6 7^ 0, which would 
yield no dark matter asymmetry, according to Eq. (13). But asymmetries 
in other species of particle can be generated in an analogous manner. For 
example, suppose that the discrete Z^ symmetry in Table I is Z$ and there 
exist complex scalar fields S which are SU (6) singlets and transform under 
Z 8 as S — > uj 4 S = —S. Then, by Table I, one sees that the S can have the 
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Yukawa coupling ((S and also have explicit superheavy masses from terms 
of the form M$S*S and AgS S + H.c. The out-of-equilibrium decays of the 
superheavy S particles can generate a ( — ( asymmetry. 

We conclude by observing that the model presented here is not unique, 
but is meant to illustrate the general point that grand unification, especially 
with large unitary groups, entails the existence of Standard-Model-singlet 
fermions that could be the dark matter. Such dark matter particles would 
be very difficult to detect. For example, in the model presented here, the 
dark matter particles have no Standard Model couplings except 0(v 2 /M 2 ) 
couplings to the Z boson through the mixing shown in Eq. (4). It is char- 
acteristic of these scenarios that there will be extra Z bosons at energies 
near the weak scale, and that the dark matter will couple to these extra Z 
bosons. Indeed, the primary means by which the dark matter particles would 
be produced in accelerators would be via the production of extra Z bosons 
and their subsequent decay into dark matter particle-antiparticle pairs. 
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